UNIQUENESS RESULTS FOR THE PHASE RETRIEVAL PROBLEM OF 
FRACTIONAL FOURIER TRANSFORMS OF VARIABLE ORDER 



PHILIPPE JAMING 

Abstract. In this paper, we investigate the uniqueness of the phase retrieval problem for 
the fractional Fourier transform (FrFT) of variable order. This problem occurs naturally in 
optics and quantum physics. More precisely, we show that if u and v are such that fractional 
Fourier transforms of order a have same modulus |J-" a it| = \T a v\ for some set r of a's, then 
v is equal to u up to a constant phase factor. The set r depends on some extra assumptions 
either on it or on both u and v. Cases considered here are u, v of compact support, pulse 
trains, Hermite functions or linear combinations of translates and dilates of Gaussians. In 
this last case, the set r may even be reduced to a single point (i.e. one fractional Fourier 
transform may suffice for uniqueness in the problem). 



1. Introduction 

Usually, when one measures a quantity, due to the nature of measurement equipment, 
noise, transmission in messy media... the phase of the quantity one wishes to know is lost. 
In mathematical terms, one wants to know a quantity ip(t) knowing only \<p(t)\ for all t E R rf . 
Stated as this, the problem has too many solutions to be useful and one tries to incorporate 
a priori knowledge on ip to decrease the under-determination. Problems of that kind are 
called Phase Retrieval Problems and arise in such diverse fields as microscopy (see e.g. \17\ 
E31 ESI HU E01 EE]), holography [2D [57], crystallography [Ml [52], neutron radiography [3|, 
optical coherence tomography [53 1, optical design [20J, radar signal processing [32J, quantum 
mechanics p~3l HH [32j [301 ESI EE] to name a few. We refer to the books [22J [55], the review 
articles [331 H31 E21 [35] for descriptions of various instances of this problems, some solutions 
to it (both theoretical and numerical) and for further references. 

The particular instance of the problem we are concerned with here deals with the so- 
called Fractional Fourier Transform (FrFT). Let us sketch a definition of this transform that 
is sufficient for the needs of the introduction (a precise definition follows in Section [3.3. 
First, for u G L 1 (M d ) n L 2 (R d ) we define the Fourier transform as 



•Fu(f )=«(£)= / u(t)e- 2iw{t & dt, 



and then extend it to L 2 (M. d ) in the usual way. Here and throughout the paper | • | and (•, •) 
are respectively the standard Euclidean norm on M. d and the corresponding scalar product. 
The inverse Fourier transform is denoted by J 7-1 . For a 6 R \ 7rZ, we define the fractional 
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Fourier transform of order a via 

T a u(0 = c a e- i7T ^ 2cota F[e-' l7T ^ 2cota u](Z/sma). 

where c a is a normalisation constant. We define .Fou(£) = u(£), J r 7r w(^) = u(— £). Also note 
that, = ^> = F an d that T a Tp = J- a +/3- 

This transform appears naturally in many instances including optics |48| . quantum me- 
chanics |38[ 144]. signal processing (5JHH]... We will detail below several instances where the 
fractional Fourier transform occurs and where one is further lead to the question of recov- 
ery of a function u G L 2 (R) from the phase- less measurements of several fractional Fourier 
transforms {l-FauHaer- More precisely, we deal with the following question: 

Problem 1 (Phase Retrieval Problem for the fractional Fourier transform). 

Let u,v G L 2 {M. d ) and let r C [0, 7r) be a set of indices (finite or not). Assume that \J r a v\ = 
\T a u\ for every a G r. 

(i) Does this imply that v = cu for some constant c G C, |c| = 1? 

(ii) If we restrict u G T> for some set T> C L 2 {W d ) do we then have v = cu for some 
constant c G C, |c| =1? 

(iii) If we further restrict both u,v G P for some set T> C L 2 {M. d ) do we then have v = cu 
for some constant c G C, \c\ = 1 ? 

In the first two cases we say that u is uniquely determined (up to constant multiples or up 
to a constant phase factor) from {l-F^-u^a G r}. In the last case we say that u is uniquely 
determined (up to a constant phase factor) from {\J- a u\,a G r} within the class T>. 

The usual phase retrieval problem is the case r = {vr/2} {i.e. when T a = J- is the usual 
Fourier transform) within the class T> of compactly supported functions or distributions. A 
part from this, the most famous problem of this sort is due to Pauli who asked whether \u\ 
and \Tu\ uniquely determines u up to constant phase factors i.e. here r = {0,7r/2}. Several 
counter-examples to this question have been constructed {see e.g p~H[30l[32] an d [13] for the 
state of the art on the problem). In order to construct those counter-examples, it was shown 
that u is not uniquely determined from \u\ and \Fu\ within the class {^finite a n( x ~ x i)} 
where 7 is the standard Gaussian nor within the class {^fi nitc o«X[o,i/2] { x ~ *)}• We will show 
that for a well chosen, \u\ and \ J- a u\ uniquely determines u up to constant multiples within 
these classes. 

Note that Reichenbach [51] conjectured that there is a unitary operator U on L 2 (IR c( ) such 
that I f\, I /| and \Uf\ uniquely determine / up to a constant phase factor. Our results thus 
show that the fractional Fourier transform is a good candidate. 

Further, we show that if the set of indices r = [0, 7r), then \J- a u\ , a G r uniquely determines 
u up to constant multiples and that the set r can be reduced to a discrete set when u and 
v are compactly supported. Moreover, we provide theoretical reconstruction formulae in this 
case. The numerical aspects of those formulae as well as further algorithms are postponed to 
forthcoming work. 

This paper is organised as follows. We start with a section in which we present various 
instances of the fractional Fourier transform in physics. The following section is devoted to 
preliminaries on Fourier analysis and complex variables. We also solve the phase retrieval 
problem for one FrFT phase-less measurement there in the class of compactly supported 
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functions. We then devote Section to two particular cases of phase-less measurements of 
windowed transforms (Fourier and wavelets). Section [5] is devoted to the solution of Problem 
[Hand is divided in several subsections dealing with the three aspects of the problem. 

2. Physical models behind the Phase Retrieval Problem for the fractional 

Fourier transform 

In this section, we will present three instances of the phase retrieval problem for the 
fractional Fourier transform. The first and main one stems from optic measurements in the 
Fresnel domain. We will here sketch how the fractional Fourier transform arises there. We 
then quickly present two more instances of our problem stemming from quantum physics. 

2.1. Diffraction phenomena in the optical far field. An excellent derivation of how 
the fractional Fourier transform (or more precisely the Fresnel transform) occurs in optics 
can be found in [35], Section 3.1]. More detailed accounts can be found in [lOj [24j [54] . We 
will here sketch the main features of the way the FrFT appears in optics, leaving aside full 
mathematical rigour for which we refer to the above mentioned texts. Roughly speaking, 
the fractional Fourier transforms are adapted to the mathematical expression of the Fresnel 
diffraction, just as the standard Fourier transform is adapted to Fraunhofer diffraction. This 
connection seems to have been made for the first time by Pallat-Finet [49] . Let us now switch 
to the physical presentation. 

Light is an electromagnetic wave with coupled electric and magnetic fields travelling 
through space. In an homogeneous isotropic medium, like free space or a lens with con- 
stant refractive index, the electric and magnetic field vectors form a right-handed orthogonal 
triad with the direction of propagation. Disregarding polarisation, the field can be described 
by a scalar function U(x;t) representing either the electric or the magnetic field amplitude. 
As the light used generally exhibits a strong monochromaticity, the time dependence of the 
field is a harmonic one and can thus be explicitly written as 

W(x;t) = R[U(x)e- lujt ]. 

Here, oj denotes the angular frequency of the light, and the complex- valued amplitude or 
"phasor" J7(x) depends on the spatial coordinates x = (x, y, z) € M 3 only. If both quantities 
U(x;t) and U(x) represent an optical wave, they must satisfy the Helmholtz equation 

(2.1) ^A- ij^Wx;i) = and (A + k 2 ) C/(x) = 0, 

£)2 £)2 pil 

where A = + + J-^ is the usual Laplacian. The Helmholtz equation directly follows 
from the Maxwell equations under the condition of a homogeneous medium and in absence of 
sources. The quantity k is termed the wavenumber or propagation constant of the medium 
and is related to the light velocity c, the angular frequency u, and the vacuum wavelength A 
by 

uj 2-nv Co 
k = — = — — and c = — , 
c A v 

whereby v is the refractive index of the medium, e.g., in vacuum Uq = 1, and cq is the light 

velocity in vacuum. 
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Figure 1. Diffraction by a planar object 



We will now study diffraction effects occurring from a planar object placed in an opening 
of a plane screen P. as illustrated in the previous picture. In order to do so, we will appeal 
to Green's theorem: 

(2.2) - J s ^(x)^(x) - I7(x)^(x) j dx = |jf L>(x)AC/(x) - C/(x)AC/(x) dx 

where S is the boundary of a domain f2 and is smooth closed and orientable and where n is 
the inner normal to that boundary. 

The integration surface S is segmented into three disjoint parts, i.e., S = A U B U C. The 
boundary A is chosen across the object location in the screen plane P, B is the opaque part 
of the screen B = P \ A, and C is the boundary of a half-sphere containing the observation 
point xo. 

Let us now consider the unit-amplitude spherical wave expanding about the observation 
point £: 

eJ*olK-x|| 

and note that 

(A + fc 2 )G (x;£) =47n%_ x . 

Let xi be the symmetric of xo with respect to P so that ||xo — x|| = ||xi — x|| for all x £ P 
and define 

f/(x) = G (x;x ) - G (x;xi). 
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Then U satisfies the following conditions 

(2.3) (A + k 2 )U(x) = 4tt(5 x _ X0 - <5 x _ Xl ) in O 

(2.4) U(x) = on P 

(2.5) ||x — xo || ~ \ —s-0 as ||x — xq|| — > oo. 

Prom this, and the fact that U also satisfies the Helmholtz equation, we derive from (|2.2p the 
so-called integral theorem of Helmholtz and Kirchhoff (see e.g. [10\ p. 377]): 

(2-6) ^) = - i// |'^M ; ?( x )_ U(x )^M) dx . 

Now, while U is identically on the plane P, its normal derivative is not. The solution of the 
diffraction problem is found by specifying Kirchhoff boundary conditions on A and B, and 
the Sommerfeld radiation condition on C: 

on A: tf(x) = I7.(x) and ^(x) = ^(x) 

on on 

dU 

on B : Uhc) = and — (x) = 

on 

on C : lim ||x|| ( -— ^(x) — ikU(x) ] = 




>+oo \ dn 

whereby U s (x) is the field incident on the screen. Let us comment on these assumptions. The 
Kirchhoff boundary conditions are an idealization to the real field distribution on the screen 
by assuming that the field and its derivative across A is exactly the same as they would be in 
the absence of the screen. In the geometrical shadow across B the field is simply set to zero. 
Although these assumptions are reasonable, they fail in the immediate neighborhood of the 
rim of the opening. The Sommerfeld radiation condition on C is satisfied, if the disturbance 
U(x) vanishes at least as fast as a spherical wave. 

Integrating these boundary conditions into the integral (|2.6p and letting the radius of the 
half-sphere C go to infinity, ()2.6p yields 

(2.7) VM—If/VM*®*. 

A 

Now, a simple computation shows that 

dU _ exp(2fc||x -x ||) / 1 \{n,x-x ) 



dn ||x — xo|| V ||x — xo||y ||x — xo 

.exp(ifc[|x - x ||) (n,x- x ) 

— 2/m rr- 



||x-x || ||x-x || 

if ||x — xo|| » A. 

At this stage, it is useful to introduce the small angle approximation wherein the angle 

between n and x — xo is small, that is ^tA [p- ~ 1. For this, we establish reference 

x - x o|| 
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coordinates (xi, X2, X3) relative to the plane P, centered on the region A. We take the X3 axis 
to be perpendicular to the plane P. Further, we assume that the measurement plane in which 
xo := (^1,^2)^3) lies is far from the screen plane P. In other words, if x = (xi,X2,x^) G P, 
then X3 = and ||(xi — £i,X2 — £2, 0)|| <C £3- But then 

(6-^l) 2 , (£2-S2) 2N V2 



X || = 6(1 + ^ + £2 



3 

2 



(d-xiY , ({2-12) 



Finally, under the various approximations made so far, 
dU exp(ik£s) f ik , 2 



exp(^|((ei-x 1 ) 2 + (6-X2) 2 )) 



dn iA6 

It follows that the field U at an observation point xq far enough from the screen is expressed 
in terms of the field U s on the screen as 

^ (xo) = -^Slr / 1 Us{x) exp Gl - xi? + (6 - X2)2) ) dx - 

A 

A simple computation then shows that 
k 

with cot a 



2vre 3 

Further occurrences of the fractional Fourier transform in optics may be found e.g. in 
[Ml EH SOI S3 HHJ [66] and in the references therein. 

As optical measurement devices are not sensitive to phase, we are lead to the question of 
reconstruction of U s from the measured quantity |f/(xo)|; or equivalently 

2AV4Al + k^ - ar§ COt k ' (2<Z ] [Us] (a ' 6) ' ' 
This leads us naturally to Problem [TJ which we rephrase as follows: can one reconstruct a 
function f from the phase-less measurements \F ai f\,... \F aN f\ with < ct% < ■ ■ ■ < ajv- 

The problem of a single measurement seems not to have attracted any specific attention. 
Indeed, the reader will easily check that this reduces to the case a\ = ir/2 i.e. the usual 
phase retrieval problem for the Fourier transform, thus a full description of the solutions can 
be given in terms of the so-called zero-flipping phenomena. For sake of completeness, we 
have included a detailed proof in Section 13.3.21 below. 

Generalisation to the case in which the intensity of an object and that of its fractional 
Fourier transform, or the intensity of any two of its fractional Fourier transforms, is known 
has been addressed in [19[ l42l [T6| [T2] to name a few. However, the problem of retrieval from 
more than two Fresnel transforms does not seem to have received anywhere near the attention 
received by the two- intensity problem and we are only aware of the papers [111 l46l El El [T], [18] . 
The problem of recovery from multiple-intensity observations related through general linear 
transformations is discussed in I3T1. 
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However, from the mathematical point of view, those papers do not present a rigorous proof 
that the object can be reconstructed from phase-less measurements, but mainly show, using 
numerical exploration, that an adaptation of the Gerchberg-Saxton algorithm will provide a 
reasonable solution. It is our aim here to provide a complete mathematical proof of facts of 
that nature. 

2.2. Quantum Tomography. Recall that the Wigner transform is defined for u G L 2 {W) 
by 

/+oo 
u{x + t/2)u(x - t/2)e 2i ^ t dt. 
-oo 

An important property of the Wigner transform is that its marginals are given by: 

x> r+oo 

W(u)(£,x)dx = \u(£)\ 2 and/ W(u)(£, x) d£ = \u(x)\ 2 . 



Because of this property, this transform was proposed as a joint time-frequency probability 
distribution with marginals the distributions of position and velocity probability. Unfortu- 
nately W(u) always takes negative values, unless u is a Gaussian, so that this interpretation 
has some defect. 

Using the link between the Wigner distribution and the ambiguity function and the fact 
that the ambiguity function of a fractional Fourier transform (see Section \3. 41 for details), one 
can check that marginals in other directions are given by the modulus of a fractional Fourier 
transform. More precisely: 

/+oo 
W(u)(t cos a — r sin a, t sin a + r cos a) dr = \F a u(t)\ 2 . 
-oo 

In particular, if we introduce the Radon transform of a function F G L 1 (IR 2 ) n L 2 (R 2 ) on M 2 
as 

lZ[F](a,t) = I F(x)dx = / F(t cos a — r sin a, t sin a + r cos a) dr 

J x :(x,(cos «,sin ot))=t Jl 

(note that this function is only defined for almost all t). Then the previous relation reads 
(2.8) K[W(u)](a,t) = \F a u(t)\ 2 . 

Because of this property, the problem of reconstructing u from \F a u(t)\ 2 for various a's is 
equivalent to the quantum state tomography problem (see [36] and references therein for 
more on this problem and |38| for the link with the fractional Fourier transform): 

Problem 2 (Quantum State Tomography). 

Let u, v € L 2 (IR) and let r C [— 7r/2,7r/2] be a set of angles. Assume that 

K[W{u)](a,t) = K[W(v)](a,t) for all a G r and all t G M. 

(i) Does this imply that v = cu for some constant c G C, |c| = 1 ? 

(ii) If we restrict u G T> for some set T> C L 2 {M) do we then have v = cu for some 
constant c G C, |c| = 1? 

(hi) If we further restrict both u, v G T> for some set T> C L 2 {U. d ) do we then have v = cu 
for some constant c G C, |c| = 1? 
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It follows that uniqueness and reconstruction results from partial Radon data may be 
transferred into similar results for the Phase retrieval problem. However, such techniques may 
not always be applicable. The main issue here is to have sufficient "concentration" of W(u) 
near the origin so that it may be reconstructed via inversion of its Radon transform. Such 
concentration is unfortunately not always available, for instance, W{u) is never compactly 
supported, unless u = 0, so that most uniqueness results for the Radon transform do not 
transfer directly into uniqueness results for the Phase Retrieval Problem for fractional Fourier 
transforms. 

A further issue is that of practical recovery by appealing to reconstruction algorithms for 
the Radon transform such as filtered back projection (see e.g. [13 Section V] or [50j Section 
3.6]). Recall that this algorithm allows to reconstruct a function Wfbp(u) from the data 
lZ[W(u)](ak, tp), ak = kir/N, tg = i/N for some N large enough. Note that, according to 
(|2.8p . this amounts for reconstruct W{u) from the data \F ak u{t?)\. This function Wfbp(u) 
is then an approximation of a smoothed version $ * W(u) (where * is the convolution over 
of W(u), up to an error term of the form 



< sup / r|J 7 [M / (u)](r cos#,rsin#)|dr 

0e\O.2n) J r>R 



6»e[0,27r) Jr>R 

where R depends on N, the "sampling scale" . Of course, there is still the issue of reconstruc- 
tion n accurately from &*W(u). Note also that a large number of "angular" measurements 
is needed here, an issue that may imply that such an algorithm may not be used in practice. 
We postpone a detailed study of this algorithm to a forthcoming paper. 

2.3. Free Shrodinger Equation. The formulation of the problem in this section corre- 
sponds to a question asked to the author by L. Vega [62]. It can also be found rather 
implicitly in |36| . 

Let us recall that the solution of the Free Shrodinger Equation 

(2.9) (id t u+±Alu = 

[u(x,0) = u (x) 
with initial data no £ L 2 (M) has solution 

u( X: t )= f e-^t+^^iO d£ = F~ x [e-^ 2 *u5] ». 
Jr 

A straightforward computation then shows that, for a S (— it/2, tt/2), 

( jJa/2 \ d 

F a uo(Z) = e~™W cota -u(C/cosQ,cota). 

y y | cos a | I 

We may thus rephrase the Phase Retrieval Problem as follows: 

Problem 3 (Phase Retrieval Problem for the Free Shrodinger Equation). 
Let uq,vq G L 2 (R d ) and let u and v be the solutions of the Free Shrodinger equation (I2.9P 
with initial value no and vq. Let r = {ti}i & i C [0, +oo) be a set of times (finite or not) at 
which one measures \u\ and \v\. Assume that 

\u(x,ti)\ = \v(x,ti)\ forall i £ I. 
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(i) Does this imply that u = cv for some constant c G C, |c| = 1 ? 

(ii) If we restrict u G T> for some set T> C L 2 (M. d ) do we then have u = cv for some 
constant c G C, |c| = 1? 

(iii) If we further restrict both u,v ET> for some set T> C L 2 (M. d ) do we then have u = cv 
for some constant c G C, |c| = 1? 

Of course, a positive answer to any of these questions implies a positive answer to the 
subsequent ones. 

Note that Pauli's problem whether an no G L 2 (U. d ) is uniquely determined (up to constant 
multiples) by |iio| and \uq\ may then be seen as a particular case of the above problem with 
r = {0,oo}. 

3. Preliminaries 

3.1. Reconstruction of an entire function of finite order from its modulus on two 
lines. In this section we gather some information about entire functions of finite order. We 
then prove that an entire function of finite type is determined up to a constant phase factor 
by its modulus on two well chosen lines through the origin. Although this result is not of 
fundamental importance in this paper and will only be applied in Section 01 it is at the 
heart of the philosophy of the paper. Indeed, the fractional Fourier transform J- a u may be 
interpreted as the restriction of a function u to the line l a = {(r cos a, r sin a) G M 2 , r G M} 
in the "time- frequency plane" (see Section I3.4p . We thus want to reconstruct u from its 
modulus on a few lines in the plane. Of course, in general {J- a u i ol G t} is not the restriction 
of a single entire function U to the lines {£ a ,a G t} so that the result presented here does 
not directly solve Problem [TJ 

3.1.1. Preliminaries on entire functions of finite order. Results in this section may be found 
in most books on one-variable complex analysis. 

Definition 3.1. Let / be an entire function (i.e. a function that is holomorphic over C). 
For r > define Mf(r) = max |/(z)|. The order p of / is defined by 

\z\=r 

log log M f (r) 
p = hm sup . 

r^+oo log r 

In other words, / is of finite order p if p is the infimum of the non-negative numbers A 
such that Mf(r) < e r e for any e > and r sufficiently large. 
Now, for k G N and £ G C, let us denote by 

1 - f\ if k = 

l-|Jexp(| + ±^--- + ±|£) otherwise 

Note that £((,e ! °z) = E(e~ ia C, z) and that E((^,z) = E(£,z). For a non-negative real 
number t we denote by [t] its integer part, that is, the integer d such that d < t < d + 1. 

Theorem 3.2 (Hadamard Factorisation Theorem). Let f be an entire function of finite 
order p, let {z n } be its zeroes, with excluded and all zeroes are repeated according to their 



E k ((,z) 
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multiplicity. Then, for every e > 0, 

( 3 - 10 ) EuT^<+°° 

n>l 

and f can be factored as 

(3.11) f(z) = z m e p ^HE [p] (z n ,z) 

n>l 

where m is an integer and P is a polynomial of degree at most p. 

Note that (|3.10p implies that the product in (|3.1ip converges uniformly over compact sets. 
Moreover, (|3.1ip shows that an entire function of finite order is essentially determined by its 
complex zeroes. 

Finally, note that in Hadamard's factorization theorem, one may replace p by any p' > p. 
3.1.2. Reconstruction of a function from the modulus on two lines. 

Theorem 3.3. Let ai,«2 6 [0)2vr[ with a\ — ot2 ^ 7rQ. Let f and g be two entire functions 
of finite order and assume that 

(3.12) \g(xe iai )\ = \f(xe iai )\ and \g(xe ia2 )\ = \f(xe ia2 )\ for all x G R 
then g = cf where c is constant of modulus 1 . 

Proof. If / has a zero at 0, then (|3.12p implies that g also has a zero at of the same 
order. Without loss of generality, we may therefore assume that / and g have no zero at 
0. For k = 1,2, let us define the four entire functions of order p, fk(z) = f(ze iak ) and 
gk{z) = g{ze lctk ) so that |<?fc(a;)| = \fk{x)\ for all This may be rewritten as 

(3.13) 9k(z)gW) = fk{z)lW) forallzeR. 

But as this is an equality between two entire functions, it is valid for all zGC. 
Now write the Hadamard factorisation of / and g as 

/( ^) = e E£W JT E [p] {z u z) J] E [p] (z h z) 

ZiEZfClZg zidZf\Z g 

and 

g( z ) = e^^ z3 n E [p] {z u z) H E [p] {z h z) 
zieZftiZg ziez g \z f 

where p is greater than the orders of / and g. Notice that (|3.13|) implies after simplification 
of the product over the common zeroes that 

fTtj***^-*)* ft E[p] (e-^z u z)E [p] (^7 u z) 
ziez f \z g 

= e ^%A^)^ J] E [p] {e-^z u z)E [p] (^7 u z). 
ziez g \z f 

From this, we get the two following consequences: 
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[P] [P] 

- first, ^ :fi l (6,<' / "' )-.•' = ^ h N (//,<'•'"' ) for k = 1,2 thus 6j = for all j / and 

i=i j=i 
K(6 ) = 9?(a ) thus 6 = a + £7, 7 G K. 

~ zi £ Zf \ Z g if and only if e 2iak ~i G Z g \ Zf for fe = 1,2 and vice versa. It results 
that (2/ \ .Zg) U (i? 9 \ i?/) is stable under the reflections with respect to the two lines e~ tai M 
and e _4a2 ]R and thus under the rotation of angle 2(a% — o^)- But this angle is not a rational 
multiple of ir thus the orbit of a point is not a discrete set. This is not compatible with 
the fact that Zf \ Z g U (Z g \ Zf) is a set of zeroes of an entire function. As a consequence 
Zj \ Z g = Z g \ Zf = that is Z g = Zf. 

Summarising, we get g(z) = e* 7 f(z) as announced. □ 

Remark 3.4. Note that it is enough to assume that \g{xe iaj )\ = \f(xe lo ' : >)\ for x G Ej where 
Ej C K is such that Eje taj is a set of uniqueness for entire functions of order p. That is, Ej 
is such that two entire functions of order p that coincide on Eje iaj coincide everywhere. For 
instance, any set of positive measure would do. 

Remark 3.5. The key point in the previous proof is the Schwarz Reflection Principle and 
the fact that the (closed) orbit of a point under certain reflections can not be included in 
a zero-set of an entire function. One may therefore extend the previous theorem to several 
variable entire functions. We refrain from doing so as we have no specific application in mind. 

Remark 3.6. In the case a\ — 02 G 7rQ, the previous proof still leads to a result if one further 
restricts the properties of the zero set of the functions / and g e.g. to be in a strip. 

Indeed, the orbit of a point in Z := Z f\Z g U (Z g \Z f) is now a regular polygon. Therefore, 
if one assumes that the zeroes of / and g are in a strip, then Z is bounded therefore finite. It 
follows that / and g differ by at most a polynomial factor / = Ph, g = Qh where h is entire 
of finite order and P, Q are polynomials. Moreover, the union of the zeros of P and Q is a 
finite union of regular polygons. 

3.2. Some facts from Fourier analysis. 

3.2.1. Fourier transforms of compactly supported functions. We will extensively make use of 
the classical Paley- Wiener Theorem: 

Theorem 3.7 (Paley- Wiener) . Let f be an entire function such that \f(z)\ < i^e 271 " 7 ' 2 ' for 
some K > and 7 > 0. // the restriction of f to the real line is in L 2 (M) then there exists a 
function u G L 2 (M) with support included in [—7,7] such that f = Tu. 

Moreover, a theorem of Titchmarsch |58j states that, 



Theorem 3.8 (Titchmarsch |58J). Let u G L 2 (]R) be a compactly supported function so that 
Tu extends to an entire function over C. Let z n = r n e ien , < r\ < r2 < . . . be the nonzero 
zeroes of Tu, counted with multiplicity and arranged according to increasing modulus. Then 



+00 j +00 ^ 

i) ^2 — = +°° but i +e < +°° f or every e > 0, 



n=l n=l 



— < +00 while y — converges (conditionally). 



71=1 71=1 
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Note that this theorem implies that Tu has infinitely many zeroes which can been seen 
directly using Paley- Wiener's Theorem and Hadamard factorization. Indeed, an expression 
like (|3.1ip can not be in L 2 (IR) if the product part is finite since P(z)e az is not in L 2 (1R) when 
P is a non-zero polynomial. 

3.2.2. The Shannon- Whittacker sampling theorem. Results in this section may be found in [7] 
and [59] and in the numerous references therein. They may also be found in numerous books 
on Fourier analysis, although the oversampling formula is sometimes slightly erroneous. 
Let us recall that the sine function is defined by 



sine x 



s§2 if x ^ 
1 if x = 



On M. d one defines sinc(x) = sinc(xi) • • • smc(xd)- 

Theorem 3.9 (Shannon- Whittacker Sampling Formula (with oversampling). Let a be a 

positive real number and let h < —. Let 7 6 C J (R ) be a non-negative even function such 

a 

that 7(£) = for |£| > 1 and normalised by 



[ 7(0d£ = l- 



Let ip G L 2 (M. d ) be such that its Fourier transform be supported in [—o~, a] d . Then ip is uniquely 
defined from its samples <p(hk/2), k G Z d by the formula 

hk\ ^ ( _ ( ^AA„-_„r,_ (_ hk\ 



(3.14) <p(x) = (2 T+ ) d ^ (y) ^ ( r - ( x " t)) sinc2?rr - 



2 7 



where r± } / 

This result is usually stated without oversampling, that is the case ah = 1 in which case 
it reduces to the usual Shannon- Whittacker Formula: 



■ 



(3.15) „(*) = (2a)- yj „ (^) sine ^ - " 

The use of oversampling improves the convergence properties of the series (|3.15p . Note also 
that various error bounds are known and can be found in text books like [251 165] : aliasing (u 
is only approximatively band-limited) jittering (the samples are only taken approximatively 
at hk), truncation (the series ()3. 15[) is truncated)... 

3.3. Preliminaries on the fractional Fourier transform. 

exp I (a — 5) 

3.3.1. Definitions. For a & K \ 7rZ, let c a = — be a square root of 1 — icota. 

V I sina| 

For u G L 1 ^) and a $ vrZ, define 

F aU (£) = c rf e -^l 2 «,t Q f u ^ e -in\t\^cota e -2in{U/^a) dt 

(3.16) = c d a e- i7r ^ 2cota T[u(t)e- lnW2cota ]^/sma) 
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while for k G Z, ^knU = u and J r (2k+i)-K u (0 = u (~0- This transformation has the following 
properties : 

(1) / FaufflTavffldt = / u(tM^d£ which allows to extend T a from L 1 (M d )nL 2 (M' i ) 

jR d J«. d 

to L 2 (R rf ) as a unitary operator on L 2 (R d ); 

(2) T a Tp = Fa+p', 

(3) if u a ^(t) = e~ 2i7T ^u(t - a), a,u G R d , then 

^>a,o(£) = F a u{Z + acosa)e- i7T ^ 2 cosas ' ma - 2inMsina ; 

and 

^a«0^(0 = F a u{^ + OJ Sin a)e in ^ 2 c°sasina+2m<u^> cosa 

Further, let us define the Hermite basis functions on R by 



Hermite functions are then defined on R rf by tensorisation 

hk u ...,k d {xi,...,x d ) = h kl (xi) ■ ■ ■ h kd (x d ). 

It is well known that (/ifc) fceN d is an orthonormal basis of L 2 (R d ) and that hk can be expressed 
as hk(x) = Hk(x)e~ w \ x \ with a polynomial of degree \k\ = k\ + • • • + k d . 
The fractional Fourier transform may alternatively be defined as 

(3-17) F a [u] = ^ kl+ ''' +kd H^h kl ,...jJh kl ,...j ln . 

(fci,...,fc d )6N d 

Remark 3.10. Note that the fractional Fourier transform T^fii] (a ^ 7rZ) is obtained from the 
ordinary Fourier transform in the following way: take u, multiply it by a chirp e _l7r '*' cota^ 
take the ordinary Fourier transform, re-scale and multiply by a chirp e' m ^ cota . As a conse- 
quence, it is often straightforward to translate results on the ordinary Fourier transform {e.g. 
Paley- Wiener and Titshmarsh's Theorems that will be used below) into results for fractional 
Fourier transforms. 

3.3.2. The Phase Retrieval Problem for the fractional Fourier transform for compactly sup- 
ported functions. In this section, for the sake of completeness, we will solve the following 
Phase Retrieval Problem: 

Problem. (FrFT Phase Retrieval Problem) 

Let a G R \ 7rZ. Given u G L 2 (R) with compact support, find all v G L 2 (R) with compact 
support such that, for every £ G R, |.F a t;(£)| = |.F a u(£)|. 

To start with, let us write u a (t) = e~ t7rt2 cota u(t) so that u a G L 2 (R) with compact support 
and T a u(0 = c a e-^ 2cota Tu a ^/ sin a). It follows that |.F Q f;(£)| = is equivalent to 

(3.18) |^a(£)| 2 = |^n Q (e)| 2 . 

Now Paley-Wiener's theorem implies that J-u a and Tv a are entire functions of exponential 
type. There are two consequences. 



11 
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— Rewriting (I3.18P as 

(3.19) Fv a (OFv a (Z) = Fu a (0^(0, 

we have an equality of entire functions, so that its validity extends from £ G R to £ G C. 
- According to Hadamard's Factorisation Theorem, Tu a may be written in the form 

(3.20) Fu a {i) = ^e a ^ +K H E 1 (C,0, 

where k u is an integer, a u ,b u G C and Z{u) := {r n e l9n ,n G N} the set of complex zeroes of 
J- a (u) (all depending on a). Moreover, Titchmarsch's theorem implies that the product in 
()3.20p converges uniformly over compact sets. A similar representation holds for J~v a . 

It follows then from ()3.19j) that k v = k u , Kea v = Rea u , Keb v = Reb u and (" G Z v if and 
only if either £ G Z u or £ G Z u . In other words, there exist uj G M, c G C with |c| = 1 and for 
each n G N, a choice e n = ±1 such that 

+oo 
n=l 

which has to be compared to (I3.20p . The convergence of this infinite product to an entire 
function of order 1 is guaranteed by Titchmarch's Theorem while the fact that one recovers 
an L 2 function is guaranteed by the fact that \J-v a \ = \J-u a \ on R. One then recovers v by 
taking the inverse Fourier transform and multiplying the result by e*^ cotQ . 

Remark 3.11. The above result is directly inspired by the result on the usual phase retrieval 
theorem (a = 7r/2) for which the corresponding result was first proved by Walter in [63], 
though in a less precise mathematical form. For proper mathematical justification of Walter's 
Theorem, we refer e.g. to |29j . 

The choice of e n = —1 amounts to changing the complex root r n e l6n of J- a u into its complex 
conjugate and is called zero-flipping in the engineering literature. 

Remark 3.12. The same proof applies as soon as we impose conditions on u and v that 
force their fractional Fourier transforms J^aftt] and ^[u] to be (extended to) entire func- 
tions of finite type. For instance, we could impose that u and v are bounded by gaussians 
|tt(x)|, \v(x)\ < Ce _A l a 'l for some A > 0. 

3.4. Preliminaries on the ambiguity function. Let us recall that the ambiguity function 
of U, V G L 2 (R) is defined by 

A(u,v)(x,y) =/,«(*+!)«(*- | )e- 2ln{t ' y) dt. 

We will simply denote A(u) = A(u,u). Note that A(u)(0, y) = F\\u\ 2 ](y). 

The reader may be more acquainted wit the following closely related transforms: the 
short-time Fourier transform, also known as the windowed Fourier transform, defined by 
S v u(x,y) = e- in ( x ^A(u,v)(x,y) and the Wi gner transform 

w( U ,v)( V ,o = j^u u+ u~ ty^^dt 
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which is the inverse Fourier transform of A[u, v) in M 2 . Further 

W(u,v)( V ,0 = 2 d A{u,v){2£,2r,) 

where v{t) = v(—t). 

Let us now list the properties that we need. They are all well-known e.g. [51 l6[ 164"]: 

(1) A(u,v) G L 2 (M. 2d ) with \\A(u, v)\\ L 2m2d^ = IM| 2 |M| 2 and is continuous, going to at 
infinity; 

(2) \\A(u)\\ oo = A(u)(0,0) = \\u\&; 

(3) A{v ) = A{u) if and only if there exists c G C with |c| = 1 such that v = cu 

(4) time-frequency shifts are transformed into time-frequency shifts: for a, u G M. d let 

u a ,u(t) = e 2iw ^^u(t - a), then 

A(u a ^v btV )(x, y) = ^{(^)Ha+b,y-u +v) ) v){x _ a + b ^ y _ UJ + ri) . 

(5) The fractional Fourier transform rotates the variables: 

A(F a u,F a v)(x,y) = A(u, v)(x cos a — y sin a, x sin a + ycosa). 
In particular, 

(3.21) A(u)(— y sin a, y cos a) = J 4(J r cr «)(0, y) = F\\F a u\ 2 ]{y). 

Property [5] was first proved in [64] but the use of the fractional Fourier transform may 
have been unnoticed as it is defined via Hermite functions. With the above expression of the 
fractional Fourier transform, this result seems to have first appeared in [44] and has been 
rediscovered several times e.g. [5| 154]. 

This property also justifies the interpretation of the fractional Fourier transform as a 
rotation of a function in the "time-frequency plane" . 

4. Reconstruction from phase-less windowed Fourier transform or wavelet 

transform with well chosen windows 

Before going on with the main scope of this paper, let us illustrate how Theorem 13.31 can be 
applied to the reconstruction of signals from phase- less measurements of particular windowed 
transforms. 

4.1. Reconstruction from a windowed Fourier transform. 

Proposition 4.1. Let j(t) = e _7r * 2 and let a G R\ irQ. Let u, v G L 2 (M.), and assume that 

(4.22) \v*l(t)\ = |w*7(*)| and |J" a [u] * j(t)\ = \F a [u] * j(t)\ 
for t £ R. Then there exists c G C with \c\ = 1 such that v = cu. 

Proof. First note that the windowed Fourier transform with window 7 is given by 

S lU (x,y) = [ u(t)j(t - x)e 2i7Tty dt 

= e -^v 2 +^y f u (t)e- n{t - z)2 dt = e- ny2+2i7TXy F u (z) 
Jm. 

where z = x + iy and F u (z) = / u(t)e~ 7T( - t ~ z ^ dt. 

it 
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Further u*7(i) = S^u(t, 0) and, using the fact that the fractional Fourier transform rotates 
the variables of the ambiguity function (Property [5] of the ambiguity function) and leaves the 
Gaussian invariant Fgj = 7 , we get that 

S jU (t cos a, t sin a) = e" i7r ' 2 cosasina A(u, 7 )(i cos a, t sin a) 
= e- i7vt2cosasina A(T a u,F a ~f)(t,0) 
= e- i7Tt2cosasina A(F a u,j){t,0) 
since J- a ^y = 7. Unwinding the computation, we obtain 

S y u(t cos a, t sin a) = e"^ 2 cosasina S 7 F a u(t, 0) = e~ int2 cosasina F a [u] * 7 (t). 

From the above computations, we see that (|4.22|) implies that |-F u (t)| = and 
\F u (te ia )\ = \F v (te ia )\ for t £ R. It remains to notice that, for u;v £ L 2 (R), F U ,F V are 
entire functions of order 2 so that Theorem 13.31 gives the result. □ 

Remark 4.2. Note that it is enough to have (|4.22p on a set of positive measure. 

Remark 4.3. By using the Paley- Wiener Theorem for Schwatz class distributions, one easily 
sees that the previous result stays true if u, v E fi'(R), provided one extends the definition of 
the fractional Fourier transform to fi' as well. 

4.2. Reconstruction from a wavelet transform. Let us recall that the wavelet transform 
with window tp £ L 2 (R) is defined for u £ L 2 (R) by 



yVV(u)(a,r) = y/a / u(t)ip(a(t - r)) dt. 
Jr 

A simple computation using Parseval's identity shows that 

W^(u)(a,T) =W^(u)(l/a,T/a) 

where 

m{u){b,0 =Vb I u(t)J(bt)e 2i ^ dt. 
Jr 

Now, let p, (5 > 0, and let us define ip(x) = ip±(x) = — so that ip+(x) = 

27r(p=fix) 

X[o,+oo) e ~ 2lTpt an d ip-(x) = X(-oo,o] e27Tpt ■ ^ ^ is supported in (—00, (3] — resp. in [— /3, +00)— 
then 

(4.23) F u , + ® := V^ + (u)(l,£) = j\(t)e- 2 ^e 2 ^ dt 

- resp. F Mj _(£) := W^ + (u)(l,£) — is an entire function of order 1. But, if a £ [0,7r/2) and 
if £ £ R such that 1 + p^ 1 ^ sin a > 0, 



F u , + (£e ia ) = / u{t)e- 2wp ( 1+p ^sinaJteMirtCcoBo^ 
Jo 

= (1 + p^S, sina) _1/2 W^(n)(l + p~ l £ sina,£cosa). 



UNIQUENESS IN PHASE RETRIVAL OF FRET'S 



17 



Now let v be such that v is also supported in (—00, 0\ and assume that |W^ + (v)(l, t)| 
|W^ + (n)(l,r)| and, if 1 + p^ 1 ^ sin a > 0, 



p p£ cos a 



p + £ sin a p + £ sin a 



/ p£cosa 



p + £ sin a p + £ sin a 



Expressing this in terms of W and then F as above, one easily checks that \F V)+ (£) | = | (£ ) | 
for £ G E and |F„ >+ (£e iQ )| = \F U:+ (£e ia )\ for £ such that 1 + p-^sina > 0. It follows from 
Theorem 13,31 and Remark 13.41 that. F v + = c+F Ut+ where |c+| = 1. But then, taking Fourier 
transforms, we obtain 

fi(*)X[o,+oo)(*) = e 2 ^[F Ui+ ](t) and 6(t)x[o,+oo)(t) = e 2 ^ l F[F Vt+ ]{t) 

thus VX[0,+oo) = C+^X[0,+oo)- 

Replacing ip + and -0- and assuming that u, v are supported in [—/3, +00) we obtain 

vX(-oc,o] = c-ux(-oo,o] with |c_| = 1. 

Summarising (and simplifying the condition), we obtain the following: 

Theorem 4.4. Let p > and Zei tb± be defined on R as i/j±(x) = — 7 c- Let u, v G L 2 (M) 

Ztawe compactly supported Fourier transforms (i.e. be band-limited) and such that \W^ ± (v)\ = 
\W^ ± (u)\. Then there exist c± G C with \c±\ = 1 suc/i f/iaf 



«(0 



c + u(£) on M+ 
c_n(£) on M~ 

Remark 4.5. Note that if we further assume that u,v are continuous, then c+ = c_. 

5. Reconstruction of a function from multiple fractional Fourier 

transform intensities 

5.1. Reconstruction from the modulus of many fractional Fourier transforms. 

We will now prove the following result: 

Theorem 5.1. In the following cases, exact reconstruction can be obtained. 

(1) Let u, v G L 2 (R) such that, for every a G [— 7r/2, it/2], \F a v\ = \F a u\. Then there 
exists c G C with \c\ = 1 such that v = cu. 

(2) Let t C [— 7r/2,7r/2] be either a set of positive measure or a set with an accumulation 
point «o 7^ 0. Let u,v G L 2 (IR) with compact support such that, for every a G r, 
\F a v\ = \F a u\. Then there exists c G C with \c\ = 1 such that v = cu. 

a 2 

(3) Let a > and define (ojfcWz by ao = ir/2 and, for k G Z \ {0}, = arctan — . 

k 

Then, given u, v G L 2 (M) with compact support included in [—a, a], if\F ak v\ = \F ak u\ 
for every k G Z, then there exists c G C with \c\ = 1 suc/i f/mi t> = cu. 

Proof. From (|3.21[) . if |.F a i;| = \F a u\ for every a G [— 7r/2, vr/2], then, from (|3.21|) . we get 
A(v) = A(u). Property [3]of the ambiguity function then implies that v = cu for some complex 
number with |c| = 1. 

Remark 5.2. As A(u) is continuous, it is enough to have a in a dense subset of [— 7r/2,7r/2] 
for the previous result to hold. 
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For the second part of the theorem, (j3.21[) implies that 

A(v)(— ysina,ycosa) = A(u)(— y sin a, y cos a) 

for every y and every a £ t. In particular, if we fix x and denote by 

t x := {—x/ tana, a G r}, 

then A(w)(x,£) = -A(u)(x,£) for every £ £ t x . But, if r is of positive measure (resp. has an 
accumulation point ao 7^ 0), then so is t x (resp. has an accumulation point at —xj tanao ~~ 
at if ao = ±ir/2) . On the other hand A(u)(x,£) is the Fourier transform of (p x defined 
by <fix{t) = u (t + § ) u (t — § ) . This function is of compact support so that A(u)(x,£) is 
an entire function in the £ variable. The same is true for A(v)(x, £). Finally, if two entire 
functions agree on a set of positive measure (resp. with an accumulation point), they agree 
everywhere. 

For the last part of the theorem, note that <p x is an L 1 function supported in the interval 
when Ixl < 2a and is for Ixl > 2a. 



\x\ \x\ 



But then, from the Shannon- Whittaker Formula, the Fourier transform A(u)(x,y) of <p x 
may be reconstructed from its samples. More precisely, 

.... . , . v-^ I h x k\ . 2tt ( h x k\ 

A(u)(x,y) = MV) = ^ [—) smc r \y + — ) 

a 1 \( h x k\ . 2ir ( h x k 



(5.24) = ^ A(u) \x, ~y j sine— I y + 

provided \h x \ < j — rr- A similar expression holds for A(v)(x, y). 

a — \x\/2 

Now recall from (|3.2ip that 

.F[|.F_ a u| 2 ](f) = A(T-av)(0,^) = A(u)(-£sma,£cosa). 
In particular, if we chose h x of the form h x = "fx then, 



2 



where = arctan — (ao = —ir/2). 

As we assumed that \F- ak v\ = \F- ak u\ for every k, it follows that A(y) {x, ^§^) = 
A(u) (x, for every x and k. From ()5.24p we then deduce that A(y) = A(u) everywhere 
thus v = cu with c G C, |c| = 1. 

It remains to choose 7 so that \yx\ < ; — r-r- for < x < 2a. For this, it is enough 

a — \x\/2 

to find the point on the hyperbola y = l/(a— x/2) at which the tangent goes through the 

2 

origin. Easy calculus shows that this point is x = a, y = 2/ a, thus any 7 < — ~ will do (see 

a 1 

Figure HJ. 

□ 
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Figure 2. The Nyquist rate and the sampling lines for A{u) 



Remark 5.3. The second part can be extended in many ways. For instance, it is enough to 
assume that u and v have fast enough decrease to force A{u) and A(v) to be entire functions 
in the second variable. For instance, this is the case if u and v have a Gaussian bound 
|«(a?)| < Ce~ x \ x \\ A > 0. 

The last part of the theorem can be improved if one increases the number of measures in 
the following way: 

Let b > a > and define (afc)fcez by ao = 7r /2 and, for k G Z \ {0}, = arctan^. 
Assume that, u G L 2 (M) has compact support included in [—a, a], and assume that we know 
\J-- ak u\ for every k G Z. Then u can be reconstructed in the following way: 

Let 7 be a function that satisfies the conditions of Theorem 13.91 that is, 7 is a smooth 
even non- negative function supported in [—1, 1] such that J 7 = 1 and let 7 be its Fourier 
Transform. Let 



ipk(x) = -Ffl-F-Q, u\ 2 } I -xsina fc + 



x\ k COS «fc 



and h x = 2x/b 2 . Then, 



b + 2x ^ , , x „ (b 2 - 2xb ( xk\ \ b 2 + 2xb f xk 



A(u)(x, y) = -^^M^l [—^~ [y-v)) Sinc -^x~ [ V ~ -W 
It is then enough to invert A(u). 
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There are several issues in the practical implementation of this formula. We need to 
estimate the error when reconstructing A(u) (aliasing, jittering and truncation) and the 
error when reconstructing u from an approximation of A(u). We postpone this study to a 
future paper. 



5.2. Reconstruction of finite combinations of Hermite functions. 

In this section, we will restrict our attention to signals u G L 2 (K d ) of the form: 

Definition 5.4. Let u G L 2 (R d ) and N a non-negative integer. We will say that u is a 
Hermite function of degree N, if u is of the form u(t) = P(f)e _7r '*' where P is a polynomial 
of degree N. 

Recall that we defined the elements of the Hermite basis /ifci,..../c d m Section 13.3.11 A 
Hermite function of degree iV may thus be expressed as 

«(*) = Yl c fc i,...,fcAi,...,fcd( t ) 

ki,...,k d eN 
ki + --- + k d <N 

with c kl ...& d G C. Further, recall that the h kl ...kd s form an orthonormal basis of L 2 (M. d ) so 
that every L 2 -function is well approximated by a Hermite function of degree N, provided N 
is big enough. 

We may now prove the following: 

Theorem 5.5. Let d > 1, M, N > be integers and let a ^ (3 £l and 7 = (3 — a. Assume 

P 

that either 7 4 Qir or 7 = — tt with p, q mutually prime integers and q > min(M, N)d. 

q 

Let u, v G L 2 (W d ) be Hermite functions of degree N and M respectively. Assume that 
= {J-'aul and that {Fpvl = \J-pu\, then there exists c G C with \c\ = 1 such that v = cu. 

The condition on 7 can also be reformulated as follows: none of e*- 77 , j = 1,..., min(M, N)d 
is real. 



Proof. Without loss of generality, we may assume that min(M, N) = N. For simplicity, let 
us give the proof only in the one-dimensional case. Write 

N M 

u(t) = e- 7Tt2 Y J C k H k (t) , v(t) = e~^Y J dkH k (t) 

k=0 k=0 

with cjv 7^ and cZm 7^ 0. It follows that, for fl£K, 

N N 

T e [u]{t) = e-^^ke~ lke H k {t) , F e [v\{t) = e^ £ e~ ik0 d k H k (t). 

k=0 k=0 
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Up to replacing Ck,dk by c^e lk/3 ,dke lk ^ we may assume that j3 = and that a satisfies 
e lja £ R for j = 1, . . . ,N. Then \v\ 2 = \u\ 2 and l-FaM) 2 = l-Fafti]] 2 is equivalent to 



(5.25) 



M 



djduH^Hkit) 

j,k=i 



N 

Y CjCkHj^Hkit) 
j,k=l 



M 



N 



djdkeW-frHjWHkit) 



j,k=l 



Looking at the highest order term in (|5.25p we obtain \du\ 2 Hf^t) on the left hand side and 
|cjv| 2 -HAr(t) 2 so that M = N and \du\ = \cn\- Up to replacing v by S^-v we may thus assume 
that djv = C N- 

Let us now look at the term of degree 2N — 1 in (|5.25p . They appear only in Hn^\Hn 
thus 



Re(d]ydN- 



Ke(cNCN-i) and Re(e ta dNdN- 



Re(e ia CNCN- 



By assumption, e ia ^ 1 so that d^d^-i = c^cn-i thus dj\r-i = cjv-i- 

Let us now assume that d^-j = c^-j for j = 0, . . . , k — 1 and let us determine djv-&. 
For this, note that the highest order term in which ctv-Zc, djv-k appear in (|5.25p is H^_kH^ 
which is of order 2N — k. As all terms of higher order have already been identified, the 
comparison of terms of order 2N — k in (|5.25|) reduces to 



Re(d N d N _ k 
By assumption, e lka ^ R so that d^d^ 



Re(c7vCAr_fc) and Re(e lfca dArci/v-fc) z 
= c N c^^k thus d N _ k 



Re(e lka c N cW^) 

' CN-k- 



□ 



Remark 5.6. One may extend this theorem so that there is no assumption on v by adding a 
third measure. More precisely, let u be a Hermite function and assume that \J- a v\ = \J- a u\, 
\Fpv\ = {Fpu] and (J 7 ^! = |.Fjjit| for three real numbers a < /3 < r] G [0, tt]. Assume that 
that < j3 — a < vr/2, tt/2 < r] — a < ir while < r/ — f3 < tt/2, then v is also a Hermite 
function. 

This is a generalised version of Hardy's Uncertainty Principle which follows immediately 
from B. Demange's proof of it [15J. The key point is that the angles a, f3, r/ define three lines 
and that none of the angular sectors that they delimit has an opening of more than tt/2. 
Similar results where considered in [27j . 

Remark 5.7. The hypothesis on 7 is necessary. For instance, one easily checks that, if G R, 
then u = hx + c^-jfiN-j and v = + CN-jh^-j satisfy \u\ = \v\ and \Tpu\ = \Fpv\. 

The proof of Theorem 15.51 may easily be adapted to prove that, under the same hypothesis, 
if \J- a v\ = \J~a u \ an d \J~pv — J- a v\ = \J-pu — F a u\, then there exists c G C with |c| = 1 such 
that v = cu. 



The second condition may be rewritten 
conjecture that the result stays 
assume that \F n v\ = \F n u\ and 



P-a 



Tnu— T a u 



P-a. 



It would be tempting to 



true if we 


let j3 — >• a 


dT g v I 




OTgU 


ae \e=a 




ae 



a in the second hypothesis, that is, if we 
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This is however not the case since, for T a u = X]^=o c j e i ^°' a j and 



06 



N 



It is then easy to check that u = + CN~ih]y-i and v = Hn + CN-ihjy-i satisfy the 
hypothesis, but v is not a multiple of u up to a constant phase factor. 
We refer to [37] for more results in this direction. 

Note that this result may also fall in the scope of Proposition 14. II provided one extends it to 
distributions since a Hermite function is the convolution of a linear combination of derivatives 
of Dirac masses at with a Gaussian. 

Remark 5.8. In [9], we defined a trivial solution of the Phase Retrieval Problem as being 
a linear or antilinear operator that sends u into a solution of the problem. In the problem 
considered here, this would be a linear or antilinear continuous operator T on L 2 (W d ) such 
that, for every u G L 2 (R d ) and every a G r, |J-* a [rit]| = |.F tt ii|. Using the density of Hermite 
functions, it is not hard to adapt the proof of (9j Proposition 3.1] to show that necessarily 
Tu = cu for some c G C, |c| = 1 as soon as r contains at least three numbers a, ft, rj satisfying 
the conditions of Remark 15.61 



5.3. Reconstruction of pulse train signals. 

ILet us now consider so called pulse train signals which are commonly used in signal process- 
ing: 

Definition 5.9. A signal u G L 2 (M) is a (rectangular) pulse-train signal (for H) if there 
exists a finite sequence (ak)kez such that 

u (t) = ^2a k X[ ,b){t - ak) = ^ a k X[ a k,ak+by 

k&L fceZ 

for some b < a/2. 

We can now prove the following: 
Theorem 5.10. Let a > and b < a/2. Let (a^), (bf.) be two finite sequences and let 



u(t) = ^2a k X[ak,ak+b) and v(t) = ^2b k X[ak,ak+b)- 

k&L fcez 



Let a G K \ S-Z. If \J- a v\ = |J r Q «|, then there exists a constant c G C with \c\ = 1 such that 



2 

V = cu. 



Remark 5.11. This result is in strong contrast with what happens when a G ^Z. In this case, 
for every k, there even exists u such that \v\ = \u\ and \J-v\ = \J~u\ has at least k solutions, 
none a constant multiple of the other. We refer to \30\ [52] for the construction. 



Remark 5.12. The result should also not be misinterpreted. Theorem 15. 101 states that a pulse- 
type signal can be reconstructed from the modulus of one fractional Fourier transform among 
pulse type signals. In general, the phase retrieval problem \T a v\ = \F a u\, u a pulse type 
signal, will have many other solutions, see Section 13.3.21 for a description of these solutions. 
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Proof. There is no loss of generality in assuming that a = 1 so that b < 1/2. An easy 
computation (see e.g. [32] ) then shows that, for x G j — -,j ' + - 

= e*™ (^a^e 2 ^) A( X[0 , 6) ) (* - j, y) 

VfcGZ / 

and a similar expression holds for A(v). According to (|3.21|) . the hypothesis of the theorem 
translates into 



a k W~je 2 ™ ky cos a j A(x[o,b) )(-ysma-j,y cos a) 



(5.26) = ^6 fe ^e 2 ^ C ° Sa U(X[o,b))(-ysina- j.ycosa) 



for all j G Z and for all y such that — ysina 6 — + *.e- for all y G I^j : = 
" 2/ • 1 2. / 1 ' 



2 sin a ' 2 sin a 
perfectly defined. 
But, for \x\ < 6, 



. Recall that we assumed that a is not a multiple of tt so that this is 



e inby 

MX[o,b))(x,y) = sinvr(6- |x|), 

L ' ' Try 

thus A(x[o,b))(~ V sin a — j, ycosa) does not vanish on a set Ij of positive measure. On Ij, 
([5.26P than reduces to 



(5.27) ^2b k b k _ je 2i7Tkycosa = 



- 2iwky cos a 



keZ k£Z 

But, as (afc),(fefc) are of finite support and as cos a ^ 0, this is an equality between two 
trigonometric polynomials. As this identity is valid on a set Ij of positive measure, it is valid 
everywhere: 

^2 b k b k ~e 2inkt = V] a k a]~e 2i7Tkt for every t G R and every j G Z. 

fcez fcez 

One then easily checks that there is c G C with \c\ = 1 such that b k = ca k for every k, thus 
u = cu. □ 

Remark 5.13. Rectangular pulse trains may be replaced by more general pulse trains u(t) = 
^Zk&z a kH{t — ak) where H is supported in [0,6) b < a/2, provided that, for every j G Z, 
A(H)(—y sin a — j,y cos a) does not vanish on a set Ij of positive measure. This is the case if 
H(t) = e^ at+ ^X[o,b)(,t), a, /3 G M since then y) = e^^^Afop^X*. y - p x /ir). 

Further, as A(H) is continuous and A(H)(0,y) = T[\H\ 2 ](y), it is enough to have 

J"[|#| 2 Kjtana) ^ for j G Z. 

This is easily seen to be the case when H = (1 — t/b)+. 
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5.4. Reconstruction of linear combinations of translates of Gaussians. In this sec- 
tion, we consider signals of the form: 

Definition 5.14. A function u G L 2 is a combination of shifted Gaussians if there exist an 
integer N, complex numbers ci, . . . , cjv and real numbers t%, . . . , tpj such that 

TV 

«(*) = 5^9j7(t- tj) 

i=i 

where 7(f) = e _7r ' 2 . 

\ 

Note that, if 5 a is the dirac mass at a, then u = I c j$tj I * 7- Thus, the next theorem 

v=i V 

has to be compared to Proposition 14.11 (extended to distributions). 

7T 

Theorem 5.15. Let a G R \ — Z and Zet u, t> 6e too combination of shifted Gaussians. If 
\J- a v\ = \ J- a u\, there exists c G C with \c\ = 1 such that v = cu. 
Proof. Let us write 

JV M 

u{t) = c jl{t ~ *j) an d v {t) = K j7(^ ~~ T j) 
i=i i=i 

with Cj,K,- G C and tj,T,- G R. Moreover the i,-'s (resp. the t,'s) are all distinct. 

An easy computation shows that A( n f)(x,y) = 2~ 1 / 2 e~ 7T ( x2+y2 ^ 2 , so that with Property 
dH) of the ambiguity function and its bilinearity, we get 

i N t \ 

A(u)(x,y) = ^J2 c j ^e- w fe +tfc ) 2 / 2 e- 7r ( ( " +ifc ^ )2+(rt ~^ )2 ) /2 . 

In particular, 

A(u)(-t sin a, t cos a) = -j= £ c .^ e ^fe+ ifc ) 2 /2 e -H(-*sin a+ifc - tj )2 + (tcos a+lifc+itj )2)/ 2 

(5.28) = ^_ ^ Ci Cfe e -fe-^) 2 /2 e -^(%e-+ tfe e— )_ 

A similar expression holds for v. 

(5.29) A( W )(-f sina,tcosa) = 2 - 1 / 2 e -' 2 / 2 £ ^ e - 7r ^-^) 2 / 2 e - l7ri (^ elQ +^ e - 8a ) . 

j,fc=l 

We will now need the two following facts: 
Fact 1. Let a G R \ — Z and let {tj}j e z and {Tj}j^z be two finite sequences of real numbers. 
If there exists j,k,j',k' G Z such that tj-e ia + t^e~ ta = Tjie lOL + Tk/e~ ta then tj = ry and 

tk = Tk'- 
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Proof of Fact 1. The condition is equivalent to 

(5.30) (tj - T f )e 2ia = r k , - t k . 

As - tk G R, (I5.30P implies that (tj - Tj')e 2ia G R. But, we assumed that a£K \ — Z, so 
that e 2ia £ R, thus t,- - = T k > -t k = 0. □ 

Fact 2. The set {e z *} z gc is linearly independent set of functions on R. 

Proof of Fact 2. Let us consider a finite linear combination of e zt that vanishes: 

N 

G(t) := A i eZj * = A j> z i G C 
3=1 

and the zj's are all distinct. Then evaluating the derivatives G,G' , . . . , G^" 1 ' at 0, we obtain 
the Vandermonde system 

Ai + A 2 + ••• + X N =0 

Xizi + X 2 z 2 + ■ ■ ■ + Xnzn = 
< . 

Mz?- 1 + X 2 z 2 N - 1 + ••• + Xnz^ 1 = 
As the Zj's are all distinct, this system has non zero determinant, thus X\ = ■ ■ ■ = Xn = 0. □ 

We can now complete the proof of the theorem. From Fact 1, each term of the form e zt 
appearing in the sum (|5.28p -resp. f|5.29p — appears exactly ones. Moreover, as \T a v\ = l^au] 
implies A(v)(— t sin a, t cos a) = A(u)(— t sin a, t cos a), Fact 2 implies that the two sums are 
equal term by term: M = N, {tj} = {tj} thus (up to renumbering) Tj = tj and then, for 
every j,k = 1, . . . , N, Kj'Kk = CjC/7. This last identity implies that there exists c G C with 
|c| = 1 such that Kj = ccj for all j and finally that v = cu. □ 

Remark 5.16. Again the result should not be over-interpreted and it only says that uniqueness 
is achieved from the measure of \J- a u\ among combinations of shifted Gaussians, provided 

7T 

a f: — Z. Again this is in strong contrast with the Pauli problem for which this result is false. 

We have only used the fact that the shifts L-'s and r,-'s are real in a mild fashion. Actually 
the same proof works if the tj 's and Tj 's are all purely imaginary or even if we restrict them 
to be in a set E £ C such that its sum-set E + E = {e + e' : e, e' G E} intersects its 2a-rotate 
only at 0: (e 2iQ (S + E)^j n (E + E) = {0}. Of course, taking tj or Tj to be complex amounts 
to taking time-frequency shifts of Gaussians. Let us now prove a more general theorem in 
which we further also allow for dilates: 

Theorem 5.17. Let a G]0, ^[ and let u,v be time-frequency translates of Gaussians, that is, 

N N d / \ M Mj 

«(*) = E E ^e 2 ^S and v(t) = £ £ ^ fc e 2 ^S 

3=1 fc=l V a j / j=l fc=l 

iy^/j Cj,fc, Kj >fc G C, ^fc, w^fc, ?7j,fc G R, and o\ > ■ ■ ■ > a N > 0, t\ > ■ ■ ■ > t m > 0. 
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Let us define the following set of angles: = ± (J^ — (xj ; and for each j for which there 
exists k such that ujj^ / (resp. rjj^ 7^ 0), let aj (resp. ^j) he defined by 

a 2 cos a + i sin a . r 2 cos a + i sin a 
e %a - ] = — -j^=^= resp. e 1 "' 3 = — 



r/ + l 



Assume that a^, a\ . . . , ajv are all distinct and that a^, 71 . . . , 7^ are all distinct. 
Then, \F a v\ = \F a u\ implies that there exists c £ C with \c\ = 1 such that v = cu. 

Proof. The proof follows in part the lines of the previous theorem. For simplicity of notation, 
let us write 



N / \ M / 

<t) = Y J C 3 e^ 1 and = £ ^e 2 ^ 7 

J = l V CT 7 / .7 = 1 V 



1/2 / x 
J7r^(w+»;)a;+(a+fe)(j/-a;+?7) J 



1/2 J 

J=l V U j / 3=1 \ T j / 

where Cj,Kj G C, crj,Tj > and Sj,tj,ojj,r]j £ R. As previously, |^"q,u| = |<F a ^| implies that 
A(v)(— t sin a,t cos a) = A(u)(— t sin a, t cos a). We will thus compute A(u) and and 
compare the expressions obtained. 

Let us first note that, writing 7< T , a ,w(i) = e 2l7TUlt j (j^^j, we obtain 

/ \ l/ 2 

A{ja,a,uj,7TAv)( X >y) = (y^J e " 

x exp ((x — a + h) 2 + ar(y — uj + f?) 2 — i(r — cr)(x — a + b)(y — u> + 77)) . 

(T + r 

In particular, J 4(7 (J)aiW ,7 Tj fe )J) )(— tsina,tcosa) is expressed in the form cexp — P a , a ,u>,T,b ,rj(t) 
where P a ,a,u,T,b,r)(t) is a polynomial of degree 2 with highest order term 

sin 2 a + ar cos 2 a — i(a — r) sin a cos a 2 

7T 1 . 

a + t 

Let us denote by C2[-X], the set of complex polynomials of degree 2. For P £ Cz[X], let 
us write J P (u) = {(j,k) : P^a^^s^ = P} and J P (v) = {(j,k) : Pr,^, % ,r fc ,t fc , % = P}. 
We may then write 

/ a . a \ 1/2 

£sina,icosa) = ^ ^ I — J Cj(^e~ p ^\ 

pgp 2 [C] (j,k)eJ P (u) + <Tk J 

A similar expression holds for A{v) with Jp(v) instead of Jp(u). We will now need an 
elaboration on Fact 2: 

Fact 3. The set of functions {e zt+< = t2 } z ^ e c is a set of linearly independent functions over R. 

We postpone the proof of this fact to the end of the proof of the theorem. As a consequence 
of this fact, we obtain that, for each P G C2[A], 

The following fact will allow us to get some information on the sets Jp{u) and Jp(v) : 
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Fact 4. Let a e]0, 7r/2[ and assume that u,v,u',v' > are such that 



sin 2 a + uv cos 2 a — z(it — v) sin a cos a 



sin 2 a + uV cos 2 a — i(u' — v') sin a cos a 



U + V 

(5.31) 

it + v 

Then v! = u and v' = v. 

Proof of Fact 4- Note that the condition on a is simply cos and sin a ^ 0. 

. . U — V U — v' l-v/u l—v'lu' 

Looking at the imaginary part m (^5.31|) gives — ■ — = — — — -, that is -n-fr- = 1 .,/,, . 

But, as one easily checks, t — > - — - is (strictly) decreasing on [0, +oo), thus one-to-one, thus 

v'/u' = v/u. 

Now, looking at the real part in (|5.3ip gives 



sin 2 a + uv cos 2 a sin 2 a + u'v' cos 2 a 



u + v u' + v' 

Factoring out u and u 1 ', we obtain 

sm 2 a _|_ / 2 a sirr^a _|_ / / / 2 a 

1 + v/u 1 + v'/u' 

As v'/u' = v/u and sina / 0, we get 1/u' = 1/u and then v' = v. □ 

As a consequence, if (j,k),(f ,k') 6 Jp{u) and (j",k") £ Jp(v), then <7j = = Tj» and 
= o"fc' = T fc" • We may therefore group all terms stemming from a given a in the expansions 
of A{u) and -A(f) and identify those terms. 
In other words, fix a > 0, and define 



j=l,...,JV 



and 



„„(«)= E <v 2 '"^(^) 



Ti=<7 



Then 
and 



j4(u)(— tsina,tcosa) = A(u CT )(— isina,tcosa) + i?„(t) 



A(v)(— i sin a, i cos a) = A(v a )(—tsma,tcosa) + i?„(i) 

where R u and -Rt, are each linearly independent both from A(u a )(— t sina, t cos a) and from 
A(v a )(— t sin a, t cos a). As A(u)(— t sin a, i cos a) = £ sin a, i cos a), this implies that 
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A(u a ){— i sin a, i cos a) = A(v a )(— i sin a, t cos a). In other words, we may now assume that 
u = u a and v = v a or, equivalently, that Uj = := a for every j,/. But then 

/— N 

(5.32) A(u)(-fsina,*cosa) = J| CjCfee - ^'^'^'"' 8 *'^^ 

i,fc=i 

and 

y Af 

(5.33) A(-u)(-tsina,icosa) = J| «; j 7^e~ F<7 ''^' 1 ' r ' t *'' 7 * ( * ) . 

We have already exploited the fact that the P<r,sj,a;j,(T,s fcl w fc 's and Pa,tj,r)j,a,t k ,Ti k ,s au have 
same second order term. Let us now use again Fact 3 (which in this case simplifies to Fact 
2) and compare those terms for which those polynomials have same first order term. 

Note that the first order term P a ,a,u),cr,b,ri(t) is 

7T (ae l{ * /2+a) + be iin/2 ~ a) - 77V 'a 4 + le ia ° + w^ 4 + le~ ia A t 

9 

o" cos cx ~\~ 1 sin 

where e* a<T = — | . But, as in the previous proof, our assumption on the angles 

implies that the two sums (|5.32|) and (|5.33p are equal term by term. In other words M = N, 
(sj,u)j) = (tj,r]j) for all j's (up to reordering) and CjCk = Kj~Kk for every j,k, thus Cj = ckj 
with |c| = 1. □ 

Before proving Fact 3, we need an intermediate elaboration on Fact 2: 

Proof of Fact 3. Let us assume that there is a vanishing linear combination with non zero 
coefficients /x^: 

(5.34) ii k e-^ kt2+Zkt = for t € M. 

k=l 

In the previous notation, we of course assume that no term is repeated z k ) / (Q, zi). As 
this is an equality of entire functions, it is true for all t £ C. 

Without loss of generality, we assume that | Ci I = max | | and that £i is real non negative 
- by changing t to te~ l( - arg ^ 1 ^ 2 . We also reorder the remaining terms to have £1 = • • • = (,m 
while all other ( k s are 7^ £1. We then re-order the M first terms to have 3ft(zi) = • • • = 
$1{zl) = max^(zfc) and K(zfe) < K(zi) for k > L. Further re-order the L first terms to have 
Im(zi) < Im (^2), . . . , Im (z£,). Note that this implies that if a > is small enough, then 
3?(zie iQ / 2 ) > ^(z k e ia / 2 ) for k = 1, . . . , M0 

Notice that £1 = max|^| implies that, for k > M + 1, 5?(Ci — Cfc) > 0- Therefore, there 
is an «o 

such that, for \a\ < a and A; > M + 1, 3?((Ci - Cfc)e iQ ) > 0. It follows that 

e -(C fe -Ci)(e^ 2 t) 2 +fe- i)eW2 t ^ Q ^ t ^ ±oo ^ 



^This means that we choose those fe's for which fa has largest modulus, rotate so that fa has also largest 
real part. Then we chose among them, those k's for which 5R(zfc) is the largest and then, after a small rotation, 
such that z\ is the unique point for which is the largest. 
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Factoring out e Ci* 2 + Z i* in (15.34j) and setting t = e* Q//2 r, we obtain 

N 

e -Cie to T 3 +a 1 e io '/ 2 r^^ e -(Cfe-Ci)e to r 2 +(2 fe -«i)e to /2 T = Q 
k=l 

Dividing by e -Ci elQ, ' r2 +^i eia/2 ' r anc | letting r — > ±00, we obtain 

M M 
(5.35) /i fe e( Zfc - 21 ) eW2T = //1 + ^ ^e( 2fe - 2l ) elQ/2T -»• when r ^ ±00 for |q| < a . 

fe=l fc=2 

As explained above, by taking a > small enoug h, M((zk - zi)e ia / 2 ) < if k = 2, . . . , M, 
thus (|5.35p takes the form \i\ + <p(r) — > where <p(r) — > as r — )■ +00. Therefore z/i = 0, a 
contradiction. □ 



6. Conclusion 

In this paper, we have studied the phase retrieval problem for single or multiple mea- 
surements of the fractional Fourier transform. This problem occurs naturally in quantum 
mechanics and in optics. 

For the single measurement case and compactly supported functions, we have shown that 
the zero-flipping phenomena occurs as for the usual Fourier transform. However, if one 
is interested in the problem for more structured signals like translates of Gaussians, then 
uniqueness is achieved. 

For multiple measurements, we have seen that uniqueness is achieved in many cases: for 
Hermite signals and pulse train signals, two measurements suffice, provided the parameters 
of the FrFT are chosen properly. For compactly supported functions, a countable set of 
measurements guaranties uniqueness and a reconstruction Formula is provided. 

In a forthcoming paper, we will explore the practical validity of this reconstruction formula. 
We will also propose a modification of the Grechter-Saxton iterative algorithm (see e.g. [22]) 
for multiple phase-less FrFT measurements and study its validity. 
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